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We derive an analytic expression for point to point correlation functions of the Polyakov loop based on the 
transfer matrix formalism. The contributions from the eigenvalues of the transfer matrix including and beyond 
the mass gap are investigated both for the 2d Ising model and in finite temperature SU (2) gauge theory. We find 
that the leading matrix element shows similar scaling properties in both models. Just above the critical point we 
obtain for 517(2) a Debye screening mass Hd/T « 4 , independent of the volume. 



1. Introduction 

The determination of the correlation length ^ 
and the screening mass iid from point-to-point 
correlation functions of the Polyakov loop is a 
non-trivial task, especially close to the critical 
point of lattice gauge theories. The difficulties 
are resulting on one hand from finite volume ef- 
fects due to the nearby transition and on the 
other hand from the unknown parametrisation of 
the heavy quark potential in the non-perturbative 
regime. 

In the transfer matrix (TM) formalism the 
levels of the transfer matrix provide an access 
to both ^ and without the introduction of 
an ansatz for the quark potential. The levels 
and matrix elements may be obtained easily from 
fits to measured plane-plane ( or zero momen- 
tum ) correlation functions since their TM form 
is known and simply exponential. In this paper 
we intend to derive the corresponding expression 
for point-to-point correlation functions. In the 2d 
Ising model we test the validity of our TM for- 
mula by comparison to the results obtained from 
plane-plane correlators. Simultaneously we are 
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able to determine, where levels beyond the mass 
gap are of importance and what can be expected 
from such an analysis. 

2. Correlation functions in the transfer 
matrix formalism 

Let us consider d— dimensional spatial lattices 
of size N'^~^L, where N denotes the number of 
points in each transverse direction and L that in 
one selected direction ( the z— direction ). The 
lattice spacing a is set to unity in the following. 
The partition funtion is then 

Z = Tt (V^) , (1) 

and V is the transfer matrix in z— direction. 
Its eigenstates | n) and eigenvalues A„ {n = 
0, 1, 2, ...) 

V|n) = A„-|n); A„ = e-^"; (2) 

^0 < Ml < < ■•■ • (3) 

are chosen to be orthonormal. In addition we nor- 
malize our partition function such that we have 
for the vacuum state 

Ao = 1, Mo = 0. (4) 
Next we define zero momentum operators by 
d{z)^N-^-J20{x^,z) , (5) 
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where 0{x±,z) is the Polyakov loop 'P{x±,z) 
for the 3 + 1 dimensional SU{2) gauge theory 
and the spin ax,z for the 2d Ising model. The 
corresponding correlation functions are 



t{z) = {6{z)-d{0)), 

resulting in 



where 



fJ-r. 



n 1 ^mn 



{n\ 0(0) |m) , 



(6) 



,(7) 



(8) 



are the level difference and the transition matrix 
element. Due to the symmetry properties of the 
eigenstates under tranformations, which change 
the sign of O, c„„ = . 

Below the critical point (3 < (3c the lowest 
nonzero energy level /^i ( the mass gap ) defines 
the large distance behaviour of the correlator. We 
therefore define the correlation length at /3 ~ /3c 

by 



(9) 



Ax (3 > (3c the mass gap /ii « and the large 
distance behaviour is given by the next level dif- 
ference A/i, so that the Debye mass is 



(10) 



where mu is the perturbative screening mass. 

A similar formula as eqj^ may now be found for 
the point-to-point correlator 



r(f) = (o(f).o(o)), 



(11) 



in the following way. The Fourier transforms of 
r(z) and T{x) are related by 



This leads us to the ansatz 



n<m 



where 
G{p;fi) 



2(1 



't^^) sinh 



4sinh2^ + Eti4sin2f 



(12) 



(13) 



(14) 



and G(p_L = 0,Pz',^J■) is just the Fourier trans- 



form of 



4L-Z) 



i.e. we have 



e + e 

added corresponding contributions for the miss- 
ing momenta components in the denominator of 
eq.^ It is now straightforward to obtain r(a;) 
by another Fourier transformation of eq.p^. 

3. Results 

We have used the 2d Ising model to test 
our ansatz, eq.^. To this end we have mea- 
sured plane-plane and point-to-point correlators 
on iV = L = 30,40,50,60 lattices. At each 
point 500000 cluster updates were performed and 
measurements taken every 10th update. In the 
twodimensional model the levels fi„ are explicitly 
known ||l|,|| . We have carried out fits to both cor- 
relators with varying numbers of levels to obtain 
the matrix elements. Both formulae lead to the 
same results, when the maximal number of levels 
is taken into account which lead to non-negative 
c^„; i.e. our ansatz is definitely confirmed. The 
final result is shown in Fig.l for N = L — 30. 
Wc find that for (3 < f3c only one term with 
Mio = /-^i, the mass gap, contributes; near (3c up 
to three terms are essential and well above (3c, 

a for 

the point-to-point correlator to SU (2) gauge the- 
ory. The Monte Carlo data ||] were computed on 
N^xA lattices, N„ = 12, 18, 26 with 10^ - 4 • 10^^ 
updates and measurements every 10th sweep. 

Here the level differences are unknown and have 
to be determined through the fit. In general we 
find a very similar behaviour as in the 2d Ising 
model. Fits with more than two levels are only 
possible on the largest lattice very close to the 
transition. Otherwise one either obtains negative 
squares of matrix elements or there is no mini- 
mum of x^. Taking into account more terms in 
eq.|l3| for the fits tends to decrease the result for 
the mass gap level. This is shown in Fig. 2. 

It is interesting to look at the behaviour of the 
next to leading level ( or level difference ) /12 ■ As 
can be seen from Fig. 3, /i2 drops from a higher 
value below (3c at the transition to a value near 
to one ( in lattice units ) and stays then rela- 
tively constant and moreover independent of the 
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Figure 1. The lowest level diflFerences ^i, ^3 and 
//21 and the respective matrix elements q , Cgp 
and c^^e"'^!^ vs. /3//3c in the 2d Ising model. 



lattice sizes used here. This second level fixes 
the large distance behaviour above /3c of the cor- 
relation functions, since /ii, as is evident from 
Fig. 2, is essentially zero there and a third level 
does not contribute. Therefore we identify it with 
IJ,£). Because we have Nr = 4 we are led to a ra- 
tio hd/T « 4, slightly higher than found with 
conventional methods Q. 

It can be shown ||] , that for iV ^ 00 

c^„ ~ TV" ; for /3 < /3e, (15) 

c^o-TV^-i; for/?>/?c, (16) 

and from finite size scaling theory for [3 (3c 

cl^ = N-'/'^-^fixN^/"). (17) 

These scaling properties are all well confirmed by 

both the 2d Ising model and the SU{2) gauge 

theory results. 
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Figure 2. Nfii resulting from one- and two- 
level fits (open and filled symbols) on = 18,26 
(squares and diamonds) lattices for SU{2). The 
inset shows for = 26 also 3-level fits (circles). 
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Figure 3. The second level 1x2 from two- 
level fits in SU{2) gauge theory for N = 
12, 18, 26 (crosses, squares,diamonds) vs. (3 . 
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